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Foreword 


This  study  is  a part  of  a seismic  research  program  conducted  by 
the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  and  spon- 
sored by  the  Project  Manager,  Remotely  Monitored  Battlefield  Surveil- 
lance System,  U.  S.  Army  Materiel  Command,  Fort  Monmouth,  New  Jersey, 
under  Project  No.  1X663719DK73  entitled  "Terrain  Penetration  and  Param- 
eter Study." 

The  study  was  under  the  general  supervision  of  Messrs.  W.  G. 
Shockley^  Chief,  MESL,  and  W.  E,  Grabau,  Chief,  En\ironmental  Systems 
Division,  MESL,  and  under  the  direct  supervision  of  Mr.  Bob  0.  Benn, 
Chief,  Environmental  Research  Branch  (ERB),  MESL.  The  computer  pro- 
grams used  for  this  report  were  prepared  by  Mr.  E.  A.  Bay lot,  Mrs.  K.  S. 
Long,  and  Mr.  J.  R.  Lundien,  ERB.  This  report  was  prepared  by 
Messrs.  Lundien  and  Benn. 

Directors  of  WES  during  this  study  and  preparation  of  the  report 
were  BG  Ernest  D.  Peixotto,  CE,  and  COL  G.  H.  Hilt,  CE.  Technical 
Director  was  Mr.  F.  R.  Brown. 
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Summary 


A study  was  conducted  to  investigate  the  generation  and  propaga- 
tion of  microseismic  signals  from  man-walking  targets.  To  illustrate 
the  relation  between  terrain  parameters  and  Rayleigh  wave  generation  and 
propagation,  the  problem  was  divided  into  four  parts,  which  were  then 
studied  graphically.  These  parts  were:  (a)  the  target-ground  interac- 

tion, (b)  energy  coupling  to  the  substrate,  (c)  Rayleigh  wave  propaga- 
tion from  the  source,  and  (d)  transmission  of  Rayleigh  waves  over  sur- 
face macrogeometry  features. 

Parameters  that  are  included  in  the  terrain  model  are:  (a)  sur- 

face rigidity  in  terms  of  nonlinear  surface  spring  constants  and 
(b)  subsurface  rigidity  in  terms  of  seismic  properties  (wave  velocities 
and  bulk  density)  and  layer  thickness.  The  terrain  combinations  include 
a wide  variation  in  site  conditions  and  have  realism  in  terms  of  environ- 
ments found  in  nature.  Sample  problems  are  included  to  demonstrate  the 
signal  construction  techniques  and  axe  analyzed  to  show  the  effect  of 
the  various  terrain  parameters  on  the  generation  and  propagation  of 
the  Rayleigh  waves. 
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GENERATION  AND  PROPAGATION  OF  MICROSEISMIC 
SIGNALS  FROM  FOOTSTEPS 


Introduction 


1.  The  vertical  component  of  a Rayleigh  wave  is  commonly  used  as 
the  energy  vehicle  to  activate  modern,  seismic,  unattended  ground  sen- 
sors (UGS).  The  amount  of  energy  generated  by  a target  of  military 
interest  (in  this  study  the  target  is  restricted  to  one  man  walking)  is 
site  dependent,  and  the  effect  of  terrain  on  Rayleigh  wave  generation 
and  propagation  is  known  to  be  substantial.  Unfortunately,  very  little 
is  known  quantitatively  about  the  relative  effects  of  specific  terrain 
parameters  on  the  phenomena.  For  this  reason,  generalizations  or  de- 
sign "theorems"  spanning  a wide  range  of  terrain  variations  cannot  be 
specified  with  confidence.  In  an  attempt  to  develop  techniques  to  study 
the  relations  between  terrain  parameters  and  Rayleigh  wave  generation 
and  propagation  the  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES)  has  divided  the  problem  of  Rayleigh  wave  generation  and  propaga- 
tion into  four  facets  that  are  amenable  to  mathematical  simulation: 

(a)  the  target-ground  interaction;  (b)  energy  coupling  to  the  substrate; 
(c)  Rayleigh  wave  propagation  from  the  source;  and  (d)  transmission  of 
Rayleigh  waves  over  surface  macrogeometry  features . 

2.  A series  of  mathematical  models  that  can  be  used  to  predict 
microseismic  signals  in  terrain  materials  have  been  developed  at  the 
WES.^  Computer  codes  have  been  written  to  implement  the  model  on  a 
Honeywell  UUO  series  digital  computer.  Predictions  can  thus  be  made 
rapidly  and  efficiently  with  the  computer  and  can  be  used  to  evaluate 
the  complex  interactions  of  the  four  facets  of  signal  generation  and 
propagation  identified  above.  The  model  can  be  used  to  study  the 
effects  of  specific  terrain  parameters  on  the  individual  facets  of  the 
Rayleigh  wave  generation  and  propagation  phenomena,  thereby  providing 
quantitative  relations  from  which  the  most  sensitive  parameters  can  be 
identified.  Further,  the  effect  of  the  sum  of  two  or  more  facets  can  be 


1 


derived.  This  flexible  tool  provides  an  effective  means  for  incorporat- 
ing consideration  of  terrain  and  its  effects  in  UGS  design  and 
deployment. 


Pui-pose  and  Scope 

3.  The  purpose  of  this  study  was  to  utilize  the  WES  computer 
model  to  study  the  effects  of  terrain,  i.e.  soil  strength,  density,  com- 
pression and  shear  wave  velocfty,  and  thickness  of  soil  layers,  on  the 
propagation  of  a Rayleigh  wave  generated  by  a man-walking  target.  The 
desired  product  was  a series  of  relations  of  terrain  parameters  and  co- 
efficients that  define  the  effect  of  each  facet  of  Rayleigh  wave  gener- 
ation and  propagation  identified  in  paragraph  1 on  the  total  analog 
wave  train.  It  was  desired  that  these  relations  be  designed  such  that  a 
graphic  solution  (i.e.  one  that  would  predict  the  amplitude  versus  fre- 
quency) could  be  generated  for  a man  walking  over  various  combinations 
of  terrain  conditions  at  any  range  from  the  source. 

k.  In  this  report,  the  theoretical  mathematical  models  used  to 
generate  the  terrain-Rayleigh  wave  propagation  relations  are  described. 
Tiie  terrain  combinations  studied  include  a wide  variation  in  site  con- 
ditions and  have  realism  in  terms  of  environments  found  in  nature.  The 
graphic  techniques  developed  are  presented  and  demonstrated  on  some 
sample  problems.  In  addition,  some  computer  predictions  are  compared 
with  graphic  solutions.  Also  included  are  several  computer-predicted 
time-domain  signals  that  correspond  to  the  computer-  and  graphic- 
predicted  seismic  signal  spectrum  curves.  Further,  some  generalizations 
concerning  the  effect  of  the  terrain  parameters  on  each  facet  studied 
are  presented. 


The  Microseismlc  Propagation  Model 

5.  The  considerations  in  the  generation  ar.vi  propagation  of  micro- 
seismic  waves  are  shown  in  fig.  1.  The  medium  in  which  the  microseismic 
waves  propagates  is  assumed  to  be  composed  of  layers  of  homogeneous. 
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viscoelastic  material.  The  source  of  the  seismic  energy  is  a force  (or 
stress)  applied  to  the  surface  of  the  medium,  which  causes  a correspond- 
ing motion  in  the  medium  as  computed  with  the  source  coupling  coeffi- 
cients. The  source  is  not  restricted  to  any  one  form  as  long  as  a 
stress  signal  is  available.  In  practical  applications  it  can  be  a 
footstep,  a vehicle  (wheeled  or  tracked),  acoustic  noise  (e.g.  aircraft 
sound  waves),  or  background  noise  (e.g.  a pinecone  dropping  to  the 
ground).  Once  the  energy  is  coupled  into  the  medium,  it  propagates 
away  from  the  source  in  various  modes . These  modes  can  be  interpreted 
as  shear  waves , compression  waves , or  the  multiple  modes  of  Rayleigh 
waves.  In  most  cases  for  intrusion  detection  devices,  the  shear,  com- 
pression, and  higher-order  Rayleigh  wave  modes  are  small  with  respect 
to  amplitude  in  comparison  with  the  fundamental  Rayleigh  wave  and  can  be 
neglected.  However,  when  these  other  modes  are  not  negligible,  their 
effect  must  be  added,  in  turn,  for  a wave  prediction.  As  the  waves 
propagate  away  from  the  source,  their  amplitudes  decrease  as  computed 
with  the  transmission  coefficients.  The  decrease  in  amplitude  is  caused 
by  geometric  attenuation  (i.e.  the  energy  is  spread  over  an  expanding 
area  as  the  signal  phase  fronts  advance  over  greater  distances)  and 
viscous  damping  (i.e.  losses  due  to  radiation  and  friction  between  soil 
particles).  If  any  topographic  feature  is  included  in  the  propagation 
path,  its  effect  can  be  estimated  with  macrogeometry  coefficients. 

Losses  from  topographic  features  are  due  to  energy  refraction,  reflec- 
tion, and  wave  conversion  (e.g.  a Rayleigh  wave  encountering  a discon- 
tinuity can  become  a source  of  new  Rayleigh  waves  and  body  waves  with 
different  transmission  coefficients).  At  the  desired  range,  the  verti- 
cal anC  radial  components  of  particle  motion  (i.e.  displacement,  veloc- 
ity, and  acceleration)  can  be  predicted. 

6.  The  prediction  equation  for  particle  motion  as  a function  of 
range  and  time  is  shown  below.  A block  diagram  of  the  calculation  pro- 
cedures is  shown  in  fig,  2. 
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where 


A = particle  motion 
l = signal  component 
p = particle  motion  number 
r = range 
t = time 

n = frequency  number 
a)  = circular  frequency 

D = Fourier  coefficients  for  source  signal 
m = mode  number 

S = surface  macrogeometry  coeificients 
B = source  coupling  coefficients 
H = transmission  coefficients 
k = wave  number 

7.  Compression  wave  velocity,  shear  wave  velocity,  mass  density, 
and  layer  thickness  sire  used  in  the  computation  of  the  wave  numbers 
and  source  coupling  coefficients  for  Rayleigh  waves.  The  computations 
are  all  made  in  the  frequency  domain,  and  the  signal  prediction  is  con- 
verted to  the  time  domain  via  an  inverse  Fourier  transform. 

8.  The  source  coupling  coefficients  and  the  transmission  coeffi- 
cients can  be  combined  to  form  the  site  transfer  function  as  shown  in 
fig.  3.  Source  coupling  coefficients  have  low  values  at  low  frequencies 
and  high  values  at  high  frequencies.  This  frequency  response  is  similar 
to  that  of  a high-pass  electronic  filter.  The  transmission  coefficients 
at  any  given  range  have  high  values  at  low  frequencies  and  low  values  at 
high  frequencies.  This  frequency  response -varies  with  range  and  is 
similar  to  that  of  a low-pass  electronic  filter.  The  site  transfer 
function  is  the  product  (or  if  the  coefficients  are  logarithmic  values, 
the  sum)  of  the  source  coupling  coefficient  and  the  transmission  coeffi- 
cient at  each  frequency.  The  site  transfer  function  has  a peaked  re- 
sponse similar  to  that  of  a band-pass  electronic  filter.  The  frequency 
at  which  the  peak  is  located  is  the  predominant  frequency  of  the  site 
and  it  varies  with  range.  Any  source  of  input  signals  with  a component 
at  or  near  this  predominant  frequency  will  be  greatly  enhanced  over  the 
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Fig.  3.  Derivation  of  the  site  transfer  function 


other  components.  This  observation  has  been  made  many  times  on  field- 
measured  data  from  various  sources . 


Graphic  Modeling  Technique 


9.  The  graphic  modeling  technique,  illustrated  in  fig.  4,  allows 
a prediction  of  predominant  frequency  of  a site  or  a prediction  of  the 
microseismic  signal  spectrum  resulting  from  a footstep  as  a function  of 
range.  Other  source  signals  can  be  substituted  for  the  footstep  signal 
if  the  stress  spectrum  is  known.  Families  of  relations  (fig.  4)  are 
derived  using  the  comprehensive  microseismic  propagation  model  dis- 
cussed in  paragraphs  5-8  over  a range  of  site  properties . The  relations 
consist  of  the  amount  jf  attenuation  of  the  source  signal  in  terms  of 
20  log  of  the  variable  in  db  ,*  as  a function  of  frequency  resulting 
from  (a)  coupling  of  target-induced  energy  to  the  ground,  (b)  propaga- 
tion of  the  signal  at  a given  range,  and  (c)  macrogeometry  effects 
(decrease  in  signal  amplitude  due  fco  refraction  and  reflection).  The 
amount  of  energy  transferred  from  the  target  to  the  ground  surface  is  a 
complex  function  of  the  target  loading  rate  and  surface  site  conditions. 
This  frequency-dependent  function  is  derived  by  a separate  computer 
program  based  on  a spring-dashpot  analog  of  the  target-ground  interac- 
tion and  is  formulated  such  that  the  nonlinear  loading  and  unloading 
of  the  ground  can  be  simulated. 

10.  The  graphic  technique  requires  the  user  to  select  a series  of 
curves  that  correspond  to  site  conditions  for  which  a prediction  is 
desired,  i.e.  a selection  of  (a)  the  source  stress  signal,  (b)  the 
site  coupling  coefficients,  (c)  the  transmission  coefficients,  and 
(d)  the  surface  macrogeometry  coefficients.  The  sum  of  these  four 
curves  (fig.  4)  (if  curves  are  of  logarithmic  values  where  sums  are 

* In  the  strictest  sense , the  term  db  is  reserved  for  the  ratio  of 
two  power  levels , Pj_  and  Pq  , as  10  log^g  P^/Pg  . However.,  for 
the  case  of  scala.r  values,  the  form  20  log^g  is  used  with  the  as- 
sumption that  the  scalar  value  has  been  normalized  with  a value  of 
unity.  Thus,  an  approximate  power  spectrum  is  available  in  db  , 
which  can  be  converted  to  absolute  scalar  magnitude  by  dividing  by  20 
and  taking  the  antilogarithm. 
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equivalent  to  products  in  the  actual  values)  produces  the  spectrum  of 
the  microseismic  signal  (particle  velocity  versus  frequency) . The  se- 
lection of  each  of  the  four  separate  curves  requires  some  prior  knowl- 
edge of  the  site  properties  or  ranges  of  values  over  which  predictions 
are  to  he  made.  Details  are  given  below. 

Source  stress  signals 

11.  At  the  present  state  of  development,  explicit  detailed  guid- 
ance cannot  be  given  for  selecting  (or  interpolating  within  the  family 
of  curves)  the  proper  stress-time  function  for  a footstep.  It  is  be- 
lieved, however,  that  the  proper  function  can  be  reasonably  inferred  if 
the  user  carefully  considers  (in  a qualitative  sense)  the  mechanism  of 
the  target-ground  interaction  in  the  light  of  how  the  curves  were 
derived. 

12,  The  source  stress  signals  generated  for  this  study  are  for 
footsteps  on  material  whose  properties  at  the  surface  may  be  different 
from  those  in  the  substrate.  Rather  distinct  superficial  stratifica- 
tion exists  for  most  natural  soils  where  the  surface  is  exposed  to  the 
environment  (rain,  freezing  temperatures,  drying  conditions,  etc.)  or 
supports  the  growth  of  vegetation.  When  stressed,  the  surface  materials 

react  nonlinearly,  even  when  subjected  to  relatively  low  stress  levels 

2 

such  as  might  be  generated  by  a man  or  a vehicle.  The  equation  used 
to  compute  these  footstep  signals  is  written  for  the  mass  of  a man  sup- 
ported by  a nonlinear  spring  (due  to  the  ground  surface)  and  dashpot 
(for  friction  and  radiation  damping) , as  follows: 

My  + C.y  + 0(y)  = F (2) 

6 o 


where 

M = effective  mass  of  man  and  soil 
e 

y = vertical  displacement 
Cg  * damping  constant 

0(y)  = spring  force  exerted  by  ground  surface 
F = applied  force  from  man's  weight 

13.  The  ground  nearly  always  shows  two  spring  forces  0(y)  : 
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that  of  the  spring  force  exerted  on  a load  under  compression  condi- 
tions, and  that  of  the  spring  force  that  acts  on  a decreasing  load 
during  rebound  conditions.  Normally  the  rebound  spring  constant  (i.e. 
shape  of  the  force-deflection  curve)  is  much  higher  than  the  compression 
spring  constant  (the  compression  spring  constant  is  often  referred  to 

as  the  coefficient  of  subgrade  reaction,  k_  , in  the  literature  dealing 

s 

o 

with  pavement  design-3) . The  compression  spring  constant  can  have  one 
of  several  shapes:  (a)  it  can  be  linear,  (b)  it  can  increase  as  the 

vertical  deflection,  y,  increases  (compaction  or  strain-hardened  condi- 
tion), or  (c)  it  can  decrease  as  the  vertical  deflection,  y .increases 
(fluidization  or  plastic  condition).  These  cases  are  illustrated  in  a 
highly  generalized  form  in  fig.  5* 


Fig.  5.  Force-deflection  characteristics  for 
nonlinear  spring  0(y) 

l4.  Forcc-deflec^ion  relations . Only  limited  information  is 
available  in  the  literature  on  force-deflection  relations  for  surficial 
materials ; nevertheless , these  data  indicate  that  the  shape  of  the 
force-deflection  curve  is  strongly  related  to  soil  strength.  Consider^ 
able  data  sure  available  on  the  strength  of  surface  soils , especially  in 
terms  of  cone  index  (which  is  used  extensively  in  vehicle  locomotion 
studies).  It  is  believed  that  the  shape  of  the  force-deflection  curve 
can  be  correctly  inferred  from  cone  index-depth  relations . To 
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illustrate,  each  class  of  force-deflection  relations,  i.e.  linear,  com- 
paction, and  plastic,  will  be  discussed  in  terms  of  cone  index  profiles 
found  in  nature. 

15.  The  linear  class  is  seldom  found  in  natural  surface  soils; 
however,  both  the  plastic  and  compaction  conditions  approach  the  linear 
condition  at  stress  levels  that  are  low  compared  with  the  shear  strength 
of  the  material.  For  the  linear  condition  to  apply  for  a man-walking 
target,  the  top  surface  layer  (i.e.  less  than  0.5  cm)  has  to  exhibit 
considerable  cohesion,  and  the  surface  cone  index  values  have  to  be  high 
(i.e.  greater  than  50).  Further,  the  cone  index  has  to  remain  high  to  a 
considerable  depth  (5  "bo  15  cm). 

16.  The  compaction  condition  is  the  most  common  one  found  in 
nature.  Typically,  soil  strength  and  density  increase  with  depth,  and, 
therefore,  a load  applied  to  the  surface  will  compact  the  surface 
against  a firmer  foundation  material.  Thus,  at  some  degree  of  compac- 
tion, the  material  can  sustain  quite  large  loads  without  further  deflec- 
tion. For  footstep  loading,  the  low  surface  strength  values  that  either 
abruptly  or  gradually  increase  with  depth  are  of  interest.  For  ex- 
ample, in  most  cultivated  fields  the  tilled  surface  soil  has  low 
strength,  but  the  plow  pan  that  often  develops  as  a result  of  the  culti- 
vation is  quite  firm.  A man's  foot  will  normally  (if  the  soil  is  moist) 
rut  the  weaker  surface  soil,  but  the  underlying  firm  soil  will  sustain 
the  man's  weight  without  deflection.  Also,  in  many  pasture  areas  where 
a good  grass  mat  and/or  root  zone  has  developed,  the  compaction  condi- 
tion will  apply.  In  this  case  the  grass  mat  will  deflect  rapidly  with 
load,  the  root  zone  will  allow  another  rate  of  deflection,  and  the 
foundation  material  will  yield  yet  another  force-deflection  relation. 

17.  The  plastic  condition  will  apply  to  those  natural  conditions 
that  permit  the  soil  material  to  "flow"  from  the  loaded  area.  Typically, 
for  man-walking  targets  these  are  areas  where  the  cone  index  values  are 
quite  low  on  the  surface  (i.e.  less  than  about  35)  and  remain  constant 

or  decrease  with  depth.  Many  areas  in  nature  can  be  found  that  exhibit 
these  characteristics,  for  example,  peat  areas  overlain  with  a thin 
veneer  of  surface  soil,  mud  flats  found  along  many  coasts  and  backswamp 


areas,  old  inland  lakes,  and  the  frost-active  zone  in  cold  regions  dur- 
ing intermittent  freeze-thaw  cycles. 

18.  Calculation  of  force-deflection  curves . Seven  cases,  se- 
lected from  previous  WE5  investigations  to  span  a wide  range  of  varia- 
tions, were  run  to  predict  the  force  applied  to  the  foundation  material 
through  the  surface  layer.  The  following  assumptions  were  made: 

&.  The  weight  of  the  man  was  160  lh  (equivalent  to  a mass 
of  72.5  kg). 

b.  The  damping  constant  changed  when  the  spring  constant  was 
changed.  The  system  was  critically  damped  under  rebound 
(tinloading)  conditions  and  0.03  critically  damped  under 
compression  (loading)  conditions. 

£.  The  unloading  spring  constant,  ks>r  , was  10  times  larger 
than  the  linear  (elastic)  approximation  of  the  loading 
spring,  kSjC  . 

d.  The  loading  rate  (i.e.  the  transfer  rate  in  one  footstep) 
was  linear  and  such  that  the  man's  mass  was  transferred 
from  one  footstep  to  the  next  in  0.01  sec. 

e_.  The  equations  for  the  spring  forces,  0(y)  , were: 

(l)  For  compaction  conditions: 


z 

max 


(3) 


(2)  For  plastic  conditions: 


0(y)  = 


ks.</ 


l+i- 

z 

max 


(h) 


where 


ks  = spring  constant  for  linear  (elastic) 
’ approximation  of  loading  spring 
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z = deflection  at  maximum  bearing  strength: 
for  compaction  condition,  zmax  is  the 
deflection  asymptote  for  y under  high 
loads;  for  plastic  condition,  zmax  is 
the  deflection  beyond  which  the  material 
appears  to  collapse. 

19.  Plates  1-7  give  (a)  the  force-deflection  curves  from  which 
the  loading  and  unloading  spring  constants  were  obtained,  (b)  the  stress 
signal  in  the  time  domain,  and  (c)  the  stress  signal  in  the  frequency 
domain.  Parts  (a)  and  (b)  are  presented  to  aid  in  selection  of  the 
proper  stress  signal  in  the  frequency  domain,  part  (c),  which  is  the 
curve  needed  for  the  prediction  of  microseismic  signals.  The  seven 
cases  are  discussed  briefly  in  the  following  paragraphs. 

20.  The  soil  conditions  for  cases  1 and  2 were  characterized  by 
high  strength  of  the  surface  soil,  i.e.  cone  index  >150.  For  case  1 
the  strength  decreased  with  depth,  whereas  for  case  2 the  strength  in- 
creased with  depth;  i.e.  cases  1 and  2 illustrate  plastic  and  compac- 
tion conditions,  respectively.  The  force-deflection  curves  for  these 

cases  are  given  in  plates  1 and  2 (left-hand  plots),  respectively.  The 

6 T 

loading  and  unloading  spring  constants  are  2.71  x 10  and  2.h2  x 10  N/m, 
respectively.  These  spring  constant  values  are  compatible  with  those 
expected  for  a well-compacted,  unimproved  road.  Comparison  of  the  two 
plates  shows  that  they  are  nearly  identical,  and  that  the  response  was 
nearly  linear,  indicating  that  both  materials  were  behaving  as  if  they 
were  on  the  low  end  of  the  plastic  and  compaction  curves . In  both  cases , 
the  applied  load  left  a rut  only  0.04  cm  deep.  As  would  be  expected 
(considering  the  force-deflection  curves),  the  stress-time  histories 
were  nearly  identical  (see  plots  in  part  b,  plates  1 and  2).  Also,  the 
force-frequency  plots  are  very  similar  throughout  the  frequency  range 
studied.  The  predominant  energy  was  generated  in  a frequency  land  from 
20  to  70  Hz. 

21.  The  force-deflection  curve  for  case  3 is  characteristic  of 
the  plastic  condition  (see  plate  3).  The  average  cone  index  to  the 
depth  of  rutting  was  in  the  neighborhood  of  25  and  remained  about  the 
same  or  decreased  below  that  depth.  The  amount  of  rutting  in  this  case 
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was  relatively  small,  i.e.  0.9  cm.  The  spring  constants  are  2.0  * 10' 
and  2.0  x 10^  N/m  for  the  loading  and  unloading  situations,  respectively. 
The  resulting  force-time  history  is  lower  in  amplitude  and  has  consider- 
ably more  energy  in  the  lower  frequencies  (less  than  30  Hz)  than  was 
found  in  the  first  two  cases. 

22.  The  spring  constants  for  case  4 (plate  4)  are  the  same  as 
those  for  case  3.  The  resulting  rut  depth  was  less  (0.67  cm  compared 
with  0.9  cm),  because  the  surface  soil  strength  was  slightly  higher  and 
increased  with  depth,  i.e.  compaction  condition.  The  force-deflection 
curve  is  almost  linear  in  the  force  range  studied,  and,  therefore,  the 
response  is  on  the  low  end  of  the  force-deflection  curve.  It  is  ex- 
pected that  this  case  can  be  found  frequently  in  untilled  areas  in 
nature.  The  stress-time  signal  for  this  case  has  a higher  amplitude 
than  for  cases  1-3;  however,  the  energy  is  contained  in  the  frequencies 
from  0 to  30  Hz  in  about  the  same  proportion  as  in  case  3. 

23.  Case  5 represents  a quite  soft  surface  soil  condition.  The 
surface  cone  index  was  about  9 and  remained  constant  or  decreased  with 
depth  (i.e.  plastic  condition).  A man  walking  over  the  surface  gene- 
rated a rut  of  9 cm.  The  spring  constants  for  this  case  are  2.0  x 10 
and  2.0  x 10^  N/m  for  the  loading  and  unloading  situations,  respectively. 
The  stress-time  signal  is  relatively  low  in  amplitude,  and  the  predomi- 
nant energy  is  contained  in  the  frequency  range  from  0 to  7 Hz  (see 
plate  5c). 

24.  The  spring  constants  for  case  6 are  the  same  as  for  case  5» 
however,  the  force-deflection  curve  indicates  rapid  initial  deflection 
under  load,  but  the  rate  of  deflection  rapidly  decreases  with  load, 
i.e.  compaction  condition  (plate  6) . This  type  of  response  would  be 
expected  if  the  surface  soil  strength  was  relatively  low  and  increased 
rapidly  with  depth.  For  this  example  the  rut  depth  was  2 cm,  and  the 
average  cone  index  to  depth  of  rutting  was  almost  20.  The  loading  and 
response  for  case  6 generated  a stress -time  signal  with  a peak  amplitude 
about  three  times  greater  than  those  of  the  first  five  cases . This 
stress-time  signal  generated  considerably  more  energy  in  the  frequency 
range  from  5 to  50  Hz  than  did  the  signal  in  the  other  cases. 
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25.  Case  7 has  the  same  spring  constants  as  cases  5 and  6;  how- 
ever, the  soil  strength  profile  increased  less  with  depth,  resulting 
in  a nearly  linear  force-deflection  curve  (plate  7).  This  loading  re- 
sulted in  a rut  depth  of  5.8  cm.  The  cone  index  was  about  15  at  the 
surface,  but  gradually  increased  with  depth.  The  stress-time  signal 
generated  by  this  target-ground  interaction  had  a peak  amplitude  of 
about  145  megadynes,  which  is  slightly  less  than  in  case  6 but  greater 
than  in  cases  1-5.  This  stress-time  signal  had  considerable  energy 
from  0 to  4 Hz.  For  frequencies  greater  than  4 Hz,  the  energy  content 
followed  the  same  trends  as  those  in  cases  3,  4,  and  5. 

26.  The  pertinent  data  for  each  case  discussed  above  are  sum- 
marized below. 


Plate 


max  Depth  Approx 


Low  plastic 
Low  compaction 
Plastic 

Low  compaction 
Plastic 
Compaction 
Low  compaction 


«/m  M/m  m cm 

< 7 

2.71 

X 

10 

2.42 

X 

10' 

0.01 

0.04 

>150 

2.71 

X 

106 

2.42 

X 

C—  ' 
O 
1— ! 

0.01 

0.04 

>150 

2.0 

X 

10  5 

2.0 

X 

106 

0.01 

0.09 

25 

2.0 

X 

10  5 

2.0 

X 

106 

0.1 

0.67 

30 

2.0 

X 

io1* 

2.0 

X 

H* 

O 

VJ1 

0.1 

9.0 

9 

2.0 

X 

10^ 

2.0 

X 

10  5 

o.o4 

2.0 

20 

2.0 

X 

io1* 

2.0 

X 

10  5 

0.2 

5.8 
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27.  As  stated  in  paragraph  11,  explicit  guidance  cannot  yet  be 
given  for  selecting  the  proper  force-frequency  curve  to  make  a predic- 
tion by  the  graphic  technique.  However,  it  is  hoped  the  discussion 
above  and  the  summary  tabulation  will  provide  insight  into  making  such 
a selection.  Ideally,  field-measured  force-deflection  curves  should  be 
used,  but  a reasonable  selection  can  be  made  on  the  basis  of  soil  type 
(soft,  medium,  or  hard)  and  footprint  depth,  e.g.  the  spring  constants 
can  be  estimated  approximately  by  noting  the  depth  of  a footprint  left 
by  a man  of  average  weight. 

Source  counlinrr  coefficients 


28.  Although  this  study  is  restricted  to  man-walking  targets, 
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the  source  coupling,  transmission,  and  macrogeometry  coefficients  can 
be  combined  to  derive  a site  transfer  function  (see  paragraph  8)  that 
will  allow  prediction  of  terrain  effects  on  a seismic  signal  generated 
by  any  source  of  relatively  low  stress  level,  e.g.  a man,  a vehicle,  a 
helicopter,  etc. 

29.  The  source  coupling  coefficient  B was  used  to  relate  ap- 
plied force  in  megadynes  to  particle  velocity  in  centimeters  per  second. 
These  relations  were  computed  using  the  WES  microseismic  propagation 
model  (equation  l)  and  are  functions  of  material  properties  (compres- 
sion wave  velocity  , shear  wave  velocity  Vg  , and  soil  mass  den- 
sity y)  and  layer  thickness  (if  layers  are  present).  Curves  for  the 
vertical  and  radial  components  of  the  source  coupling  coefficients  were 
generated  for  (a)  a homogeneous  half-space,  (b)  a layer  on  top  of  a 
homogeneous  half-space,  and  (c)  two  layers  on  top  of  a homogeneous 
half-space. 

30.  The  seismic  properties  of  materials  used  to  calculate  the 
source  coupling  coefficients  are  given  in  table  1.  They  encompass  a 
wide  range  of  variation.  For  example,  the  compression  wave  velocities 
span  a range  from  100  to  2000  m/sec.  This  range  covers  many  terrain 

5 

materials  found  in  nature,  as  evidenced  by  the  following  tabulation. 


Estimated  Compression 
Wave  Velocity 

Soil  or  Rock  Types m/sec 


Dry,  loose  topsoils  and  silts 

Dry  sands,  loams,  and  slightly  sandy  or  grav- 
elly soft  clays 

Dry  gravels ; moist  sandy  and  gravelly  soils 

Dry,  heavy  gravelly  clays;  moist  heavy  clays; 
cobbly  materials  with  considerable  sands  and 
fines ; soft  shales ; soft  or  weak  sandstones 

Water;  saturated  silts  or  clays;  wet  gravels 

Compacted  moist  clays ; saturated  sands  and 
gravels;  soils  below  the  water  table;  dry, 
medium  shales ; moderately  soft  sandstones ; 
weathered  moist  shales  and  schists 

(Continued) 


180 

to 

365 

300 

to 
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to 
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910 

to 
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i46o 

to 

1524 

l46o 

to 

1829 

IT 


Estimated  Compression 
Wave  Velocity 
m/sec 

Hardpan ; cemented  gravels  *,  hard  clays ; boulder 
tills ; compact  cobbly  and  bouldery  materials ; 
medium  to  moderately  hard  shales  and  sand- 
stones ; partially  decomposed  granites ; 

jointed  and  fractured  hard  rocks  1676  to  2438 

Hard  shales  and  sandstones ; interbedded  shales 
and  sandstones ; slightly  fractured  limestones 
and  crystalline  rocks  2438  to  3657 

Unweathered  limestones , granites,  gneiss,  and 

other  dense  rocks  3657  to  6100 

31.  The  Rayleigh  and  shear  wave  velocities  used  are  commensurate 

with  the  compression  wave  velocity.  In  many  terrain  materials  Rayleigh 

wave  velocity  can  be  assumed  to  be  equiva] ent  to  shear  wave  velocity 

and  to  be  40-50  percent  of  compression  wave  velocity.  The  mass  density 

also  spans  the  natural  range  of  variation,  i.e.  values  of  mass  density 
3 

of  1. 3-2.0  g/cm  cover  most  soil  materials  and  values  of  2.3  and  2.6 
are  typical  of  highly  consolidated  or  rock  material.  Material  G 
(table  l)  was  selected  to  represent  a frozen  medium. 

32.  The  tabulation  in  paragraph  30  provides  crude  estimates  of 
compression  wave  velocity  as  a function  of  material  type.  Compression 
wave  velocities  vary  dramatically  as  a function  of  material  type,  and 
a closer  approximation  (for  the  selection  of  the  proper  coupling  coef- 
ficient for  a specific  site)  can  be  made  if  the  empirical  equations 
described  in  the  following  paragraphs  are  used. 

33.  Prediction  equation  for  compression  wave  velocity.  A multi- 
ple regression  computer  program  was  used  to  develop  the  following  em- 
pirical equation  for  the  prediction  of  compression  wave  velocity  in 
soil  materials  as  a function  of  basic  soil  parameters.* 


* L.  E.  Link,  "A  System  for  Predicting  Seismic  Sensor  Performance,” 
Technical  Report  (in  preparation),  U.  S.  Army  Engineer  Waterways 
Experiment  Station,  CE,  Vicksburg,  Miss. 


Soil  or  Rock  Types 
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V . 
pi 


552.0  + 1* . 9H°  * 5 


W.yd. 

i_Li + 

(l-Yd./2.65) 


-79.6 


[Wj.yd.d^)] 


(5) 


where 

V . *=  compression  wave  velocity  for  i layer,  m/sec 
P1 

H.  = depth  to  bottom  of  i layer  from  surface,  cm 

* “t  h 

W.  = average  soil  moisture  content  for  i layer,  percent 

yd^  = average  dry  density  for  i^*1  layer,  g/cm^ 

e = base  of  natural  logarithm 

This  equation  was  developed  using  selected  field  data  collected  in 
Fanaru  and  West  Germany,  and  at  Fort  Huachuca  and  Yuma,  Ariz.,  and 
Vicksburg,  Miss.  It  provides  a method  for  obtaining  compression  wave 
velocity  whenever  on-site  wave  velocity  cannot  be  measured  or  estimates 
are  desired  for  seasonal  variations. 

34.  Surface  soils  are  sometimes  underlain  by  rock  instead  of  a 

second  soil  layer,  especially  in  mountainous  regions;  thus,  an  ability 
to  predict  compression  wave  velocity  for  rock  materials  is  necessary. 
However,  since  compression  wave  velocity  varies  dramatically  as  a func- 
tion of  rock  type  and  condition  (weathering,  fracturing,  etc.),  no 
attempt  was  made  to  develop  a specific  prediction  capability  for  rock. 
Instead,  it  is  recommended  that  the  last  three  materials  listed  in  the 
tabulation  in  paragraph  30  be  vised  for  estimating  values  on  the 

basis  of  rock  type  and  condition.  Although  the  tabulation  allows  only 

crude  estimates  of  V to  be  made  for  rock  materials,  it  is  believed 

P 

that  these  values  will  be  adequate  to  provide  a reasonable  picture  of 

the  substrate  because  of  the  usually  great  contrast  that  will  exist 

between  V values  for  surface  soils  and  those  for  most  rock  materials. 
P 

35.  Prediction  of  shear  wave  velocity.  A prediction  equation 

for  shear  wave  velocity  V similar  to  that  for  predicting  the  com- 

s 

pression  wave  velocity  V of  any  soil  with  depth  could  not  be 

It 
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developed  because  of  a lack  of  data.  However,  if  Poisson's  ratio  a 
is  known,  the  ratio  Vs/Vp  can  be  determined  with  the  following 
equation: 


Previous  experiments  have  shown  that  Poisson's  ratio  for  natural  soils 
normally  ranges  between  0.30  and  0.49.  By  assuming  that,  in  general, 
most  soils  can  be  described  by  a Poisson's  ratio  of  0.4,  a value  of  0.4 
can  be  computed  for  Vs/V^  or 

V = 0.4V  (7) 

s p 

By  using  this  equation  and  the  empirical  equation  for  predicting  , 
a technique  is  available  for  estimating  shear  wave  velocity  for  any 
soil  layer  with  depth. 

36.  In  rock  materials  Poisson's  ratio  may  vary  from  approxi- 
mately 0.10  to  0.40  as  a function  of  the  rock  type  heterogeneities  and 
rock  weathering.  If  an  average  value  of  0.25  is  assumed,  the  ratio 

V /V  then  becomes  O.58,  and  the  expression  can  be  written 
s p 

V = O.58V  (for  hard  rock)  (8) 

s p 

Since  all  types  of  rock  and  rock  in  every  state  of  deterioration  exist 

in  the  real  world,  it  cannot  be  assumed  that  the  expression  V = O.58V 

s p 

will  hold  true  in  all  cases.  It  will,  however,  provide  a means  for 

estimating  V for  rock  from  the  estimates  of  V shown  in  para- 
s p 

graph  30 . 

37*  Conditions  studied.  Various  combinations  of  the  materials 
described  in  table  1 were  selected  for  study.  As  stated  in  para- 
graph 29,  three  types  of  soil  models  were  used:  (a)  a homogeneous 

half-space,  (b)  one  layer  over  a homogeneous  half-space,  and  (c)  two 
layers  over  a homogeneous  half-space.  Various  layer  thicknesses  were 
used:  0.25,  0.75,  1*5,  3.0,  and  6 m.  The  major  portion  of  the  study 
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was  conducted  on  "normal"  soil  models,  i.e.  the  density  and  rigidity 
of  the  material  increased  with  depth.  However,  to  determine  the  ef- 
fects of  a frozen  layer  over  thawed  material,  a combination  of  mate- 
rial G over  material  B (see  table  l)  was  studied.  The  various  cases 
studied  and  the  plates  in  which  the  results  are  plotted  are  tabulated 
below. 


Site  Composition Source  Coupling  Coefficients 


Case 

Top  Layer 
Material 

Base 

Material 

Radial  Component 
Plate  No. 

Vertical  Component 
Plate  No. 

1 

A 

Homogeneous 

8 

9 

2 

B 

Homogeneous 

8 

9 

3 

C 

Homogeneous 

8 

9 

It 

D 

Homogeneous 

8 

9 

5 

E 

Homogeneous 

8 

9 

6 

F 

Homogeneous 

8 

9 

7 

G 

Homogeneous 

8 

9 

8 

A 

D 

(l  layer  over 
base  material) 

10 

19 

9 

A 

E 

11 

20 

10 

A 

F 

12 

21 

11 

B 

D 

13 

22 

12 

B 

E 

lit 

23 

13 

B 

F 

15 

2lt 

lit 

C 

E 

16 

25 

15 

C 

F 

17 

26 

16 

G 

B 

18 

27 

17 

B 

C-F 

(2  layers  over 
base  material) 

28 

28 

18 

B 

D-E 

29 

29 

38.  As  seen  from  the  tabulation  above,  plates  8 and  9 show  the 
radial  and  vertical  source  coupling  coefficients,  respectively,  for 
each  material  taken  as  a homogeneous  half-3pace  (cases  1-7).  The  com- 
bination shown  in  plates  18  and  27  of  a high-velocity  layer  (material  G) 
over  a base  (material  B)  is  given  as  an  example  of  a frozen  site 
(case  16).  The  rest  of  the  profiles  are  examples  of  unfrozen  sites. 

Two  examples  (cases  17  and  18)  are  given  for  a site  with  two  layers  on 
top  of  a homogeneous  half-space  (see  plates  28  and  29).  Cases  8-l6 
are  for  si+es  that  can  be  characterized  by  one  layer  over  a half-space. 
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For  many  terrain  situations  this  more  simplified  terrain  model  (i.e. 
one  layer  on  top  of  a half-space)  has  been  found  adequate  for  calcula- 
tions made  over  the  frequency  range  used  by  the  seismic  intrusion 
sensors . 

39*  A comprehensive  analysis  of  plates  9-29  to  summarize  trends 
as  functions  of  the  various  combinations  of  materials  has  not  been  made; 
however,  some  tentative  generalizations  are  presented  in  the  following 
paragraphs . 

Source  coupling 
coefficients-terrain  relations 

40.  Homogeneous  half-space.  The  source  coupling  coefficients  B 
for  a homogeneous  half-space  increase  with  frequency  and  are  propor- 
tional to  the  square  of  frequency  (see  plates  8 and  9)*  The  relative 
amplitude  of  the  source  coupling  coefficients  is  set  by  the  rigidity  of 
the  material.  For  example,  a rigid  material  is  characterized  by  a high 
shear  modulus  (a  high  shear  wave  velocity)  and  a low  tendency  to  deform 
under  an  applied  stress;  a rigid  material  has  small  source  coupling 
coefficients.  Conversely,  a less  rigid  material  will  deform  to  a 
greater  extent  under  an  applied  stress  and  has  larger  source  coupling 
coefficients.  Over  the  range  of  variation  presented  in  this  report, 
the  source  coupling  coefficients  are  approximately  proportional  to 

O 

(l/V  ) . For  homogeneous  materials,  the  radial  component  of  20  log  B 
s 

is  approximately  6 db  less  than  the  vertical  component  of  20  log  B . 
(Compare  plates  8 and  9»)  This  difference  of  6 db  means  that  the 
radial  component  of  the  source  coupling  coefficient  is  about  50  percent 
of  the  vertical  component. 

41.  One  layer  over  homogeneous  half-space  (normal  soils).  The 
value  of  20  log  B for  one  layer  over  a homogeneous  half-space  appears 
to  exhibit  asymptotic  tendencies,  i.e.  the  high-frequency  asymptote  has 
values  approximately  equal  to  the  source  coupling  coefficients  for  a 
homogeneous  half-space  composed  of  the  surface  material,  and  the  low- 
frequency  asymptote  has  values  approximately  equal  to  the  source  coup- 
ling coefficients  for  the  foundation  material  alone.  The  effect  of  layer 
thickness  on  the  value  of  20  log  B is  readily  seen  in  plates  10-27 ; in 
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general,  as  the  layer  thickness  decreases  the  coupling  coefficients  de- 
crease. Also,  as  the  layer  thickness  decreases,  the  transition  band 
bet./een  the  high-  and  low-frequency  asymptotes  shifts  to  higher 
frequency . 

42.  One  layer  over  homogeneous  half-space  (frozen  soil).  The 
curves  for  20  log  B for  a more  rigid  material  (e.g.  frozen  layer)  over 
a less  rigid  infinite  half-space  have  an  appearance  much  different  from 
that  for  the  unfrozen  case.  This  appearance  is  primarily  due  to  the 
position  of  the  high-  and  low-frequency  asymptotes,  which  are  reversed 
from  those  of  the  unfrozen  soil,  as  would  be  expected.  The  more  rigid 
surface  seems  to  be  acting  as  a plate  on  a compressible  half-space  to 
give  low  coupling  coefficient  values. 

43.  Two  layers  over  homogeneous  half-space.  The  curves  for  20 
log  B for  two  layers  over  a homogeneous  half-space  (plates  28  and  29) 
have  high-  und  low-frequency  asymptotes,  as  in  the  one-layer  case.  In 
both  plates,  the  radial  component  is  larger  than  the  vertical  component 
over  a considerable  frequency  range  and  is  characterized  by  a large 
resonance  peak  in  the  middle  of  the  transition  frequency  band.  The  ef- 
fect of  differences  in  foundation  material  shows  up  in  the  values  com- 
puted for  20  log  B . For  example,  the  situation  illustrated  in  plate  28 
is  one  that  has  4.5  m of  relatively  low  shear  wave  velocity  (i.e.  120 
and  200  m/sec,  respectively,  for  each  layer)  material  over  a highly 
contrasting  base  material  (shear  wave  velocity  = 1000  m/sec).  In  the 
case  illustrated  in  plate  29,  there  Eire  2.25  m of  low  shear  wave  veloc- 
ity (120  and  200  m/sec,  respectively)  material  over  a base  material  with 
a shear  wave  velocity  of  500  m/sec.  This  combination  results  in  the 
vertical  component  of  20  log  B being  less  for  the  situation  shown  in 
plate  28  than  for  the  situation  in  plate  29  at  low  frequencies  (less 
than  about  l4  Hz),  and  greater  or  approximately  the  same  at  higher  fre- 
quencies. The  radial  component  has  grossly  the  same  values  in  both 
situations,  but  the  rapid  change  with  frequencies  apparent  in  the  curves 
occurs  at  lower  frequencies  for  the  more  highly  contrasting  materials. 
Transmission  coefficients 

44.  The  transmission  coefficients  H are  plotted  in  plates  30-32 
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in  terms  of  20  log  H versus  tcR  , where 


and 


2rtfR 
" VR 


(9) 


k - wave  number,  m ^ 

R = range  from  source,  m 
f = frequency,  Hz 

VR  = Rayleigh  wave  phase  velocity,  m/sec 
The  coefficients  ware  generated  under  the  assumption  that  the  ground 
is  underdamped  with  a value  equal  to  0.03  of  critical  damping  (a  value 
which  is  an  acceptable  compromise  for  many  soil  conditions).  Plates 
30-32  can  be  used  to  determine  relations  of  the  site  transmission  coef- 
ficients (in  terms  of  20  log  H)  versus  frequency  provided  that  the 
relations  of  Rayleigh  wave  phase  velocity  versus  frequency  (dispersion) 
for  the  site  are  known.  V_-frequency  relations  for  specific  site  condi- 
tions  can  be  obtained  from  plates  33-42,  and  these  relations  can  be 
used  to  calculate  values  of  icR  over  the  frequencies  of  interest. 
Corresponding  values  of  20  log  H can  be  read  directly  from  plates 


30-32. 

45.  Fnase  velocity  was  computed  as  a function  of  frequency  for 
the  same  cases  as  were  the  source  coupling  coefficients  (paragraph  37) • 
The  resulting  relations  are  shown  in  table  1 and  plates  33-42  as  sum- 
marized below. 


Case 

Site  Composition 
T9P  Layer  Base 

Material  Material 

Rayleigh  Wave  Phase 
Velocity  Versus  Frequency 
Plate  No. 

1 

A 

Homogeneous 

« 

2 

B 

Homogeneous 

* 

C 

Homogeneous 

« 

4 

D 

Homogeneous 

* 

5 

E 

Homogeneous 

« 

6' 

F 

Homogeneous 

* 

7 

G 

Homogeneous 

* 

(Continued) 

* Rayleigh  wave  phase  velocity  values  are  shown  in  table  1. 
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Case 

Site  Composition 
Top  Layer  Base 

Material  Material 

Rayleigh  Wave  Phase 
Velocity  Versus  Frequency 
Plate  No. 

8 

A 

D 

(l  layer  over 
base  material) 

33 

9 

A 

£ 

34 

10 

A 

F 

35 

11 

B 

D 

36 

12 

B 

E 

37 

13 

B 

F 

38 

l4 

C 

E 

39 

15 

C 

F 

40 

16 

G 

B 

4l 

17 

B 

C-F 

(2  layers  over 
base  material) 

42 

18 

B 

D-E 

42 

46.  It  should  be  noted  that  for  cases  1-7,  i .e.  homogeneous  half- 
space, the  values  for  Rayleigh  wave  phase  velocity  do  not  change  with 
frequency,  and  the  values  listed  in  table  1 can  be  used  directly  to 
determine  the  transmission  coefficients  for  a given  set  of  site 
conditions . 

4?.  As  with  the  source  coupling  coefficients-terrain  relations, 
a comprehensive  analysis  of  plates  33-42  to  summarize  trends  as  a func- 
tion of  the  various  combinations  of  materials  has  not  been  made.  How- 
ever, layer  thickness,  as  well  as  material  properties,  has  a profound 
effect  on  the  Rayleigh  wave  dispersion  curve.  Some  of  the  tentative 
generalizations  are  presented  in  the  following  paragraphs . 

Rayleigh  wave  phase 
velocity-terrain  relations 

48.  Homogeneous  half-space.  The  Rayleigh  wave  phase  velocity  VR 

for  a homogeneous  haif-space  is  fixed  by  the  media  seismic  properties 

(i.e,,  V , V , and  y)  and  does  not  change  with  frequency  (see 
s p 

table  l).  The  ratio  VD/V  is  generally  in  the  order  of  0.90  to  0.95 

Jt\  s 

for  most  terrain  materials. 

49 • One  layer  over  homogeneous  half-space  (normal  soils).  ( See 
plates  33- 40 . ) The  Rayleigh  wave  phase  velocity  for  one  layer  over 
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a homogeneous  half-space  changes  with  frequency.  In  the  case  of  normal 
or  unfrozen  soils,  i.e.  where  rigidity  increases  with  depth,  the  value 
of  the  Rayleigh  wave  velocity  is  dependent  on  the  material  in  which  the 
majority  of  the  wave  is  propagating.  For  example , if  the  wavelengths 
are  long  compared  with  the  first-lay- *r  depth,  the  Rayleigh  wave  velocity 
approaches  the  Rayleigh  wave  phase  velocity  of  the  foundation  material. 

As  shown  in  plates  33-1*0  the  contrast  in  the  propagation  velocities 
between  the  first  layer  and  the  half-space  is  reflected  in  the  dispersion 
curves  because  the  Rayleigh  wave  phase  velocity  must  take  on  values  that 
lie  between  these  two  values.  Further,  it  can  be  seen  that  the  relation 
between  Rayleigh  phase  velocity  and  frequency  is  also  a function  of 
first-layer  depth.  In  general,  the  high-frequency  asymptote  has  values 
approximately  equal  to  the  Rayleigh  wave  phase  velocity  of  a homogeneous 
half-space  composed  of  the  surface  material,  and  the  low-frequency  asymp- 
tgte  has  values  approximately  equal  to  the  Rayleigh  wave  phase  velocity 
of  a homogeneous  half-space  composed  of  only  the  foundation  material. 

The  change  from  low  to  high  velocity  starts  at  a frequency  f where 
the  wavelength  X(X  = V_/f)  is  between  1.0  and  1.5  times  the  layer  thick- 
ness.  The  width  of  the  transition  frequency  band  (i.e.  the  width  of  the 
band  of  frequencies  over  which  the  change  occurs  between  asymptotes) 
measured  in  octaves  appears  to  increase  in  a linear  fashion  with  the  top 
layer  shear  wave  velocity.  For  example,  the  width  of  the  transition 
frequency  band  is  approximately  1-1/2  octaves  for  a top  layer  shear  wave 
velocity  of  1*0  m/sec  (see  plates  33-35)  and  3 octaves  for  a shear  wave 
velocity  of  200  m/sec  (see  plates  3 9 and  1*0 ). 

50.  Some  preliminary  comparisons  were  made  of  the  Rayleigh  wave 
phase  velocity  relations  with  the  relations  of  source  coupling  coeffi- 
cients in  an  attempt  to  derive  a means  for  collapsing  the  data  to  a 
simpler  form.  A gross  relation  for  the  vertical  component  evolved,  i.e. 

B = f 2/vjj  (see  paragraph  1*0)*,  however,  this  equation  does  not  ade- 
quately define  the  relation  between  the  two  functions,  as  evidenced  by 
the  fact  that  the  effects  of  layering  only  approximately  occur  at  the 
same  frequency.  For  example,  the  Rayleigh  wave  phase  velocity  converges 
to  form  one  curve  at  low  frequencies \ whereas  most  of  the  20  log  B 
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curves  are  shown  to  be  separated  at  these  frequencies.  This  inconsis- 
tency in  the  two  relations  is  probably  due  to  resonance  effects  (such 
as  those  caused  by  the  surface  layer  acting  on  a half-space),  and  fur- 
ther study  is  needed  before  a general  simple  equation  can  be  derived 
that  defines  the  relations  between  the  two  functions. 

51 • One  layer  over  homogeneous  half-space  (frozen  soils).  The 
Rayleigh  wave  phase  velocity  for  one  layer  over  a homogeneous  half-space 
(plate  111)  also  changes  with  frequency  for  frozen  soils  (simulated  by 
assuming  a high-velocity  layer  over  a lower -velocity  foundation  mate- 
rial) but  in  a manner  just  reversed  from  that  of  the  unfrozen  soils. 

The  same  generalizations  appear  to  hold  true  for  frozen  soils  as  for 
unfrozen  soils,  if  it  is  realized  that  the  high-frequency  asymptote  will 
have  values  approximating  the  high  Rayleigh  wave  phase  velocity  of  the 
surf ace  material,  and  the  low-frequency  asymptote  will  have  values  ap- 
proximating the  low  Rayleigh  wave  velocity  of  the  foundation  material. 

52.  Two  layers  over  homogeneous  half-space.  The  Rayleigh  wave 
phase  velocity  for  two  layers  over  a homogeneous  half-space  varies  with 
frequency  in  a manner  very  similar  to  that  for  one  layer  over  a homoge- 
neous half-space.  The  high-frequency  asymptote  has  values  approximating 
the  Rayleigh  wave  phase  velocity  of  a homogeneous  half-space  composed 

of  the  top  layer  material,  and  the  low-frequency  asymptote  has  values 
approximating  the  Rayleigh  wave  phase  velocity  of  a homogeneous  half- 
space composed  of  the  foundation  material.  The  effect  of  the  second 
layer  is  to  cause  the  transition  frequency  band  to  increase  over  that 
shown  for  a one-layer  system.  If  the  second  layer  is  thick  enough,  a 
second  asymptote  may  appear  midway  between  the  high  and  low  asymptotes, 
which  is  related  to  the  second-layer  shear  wave  velocity. 

Surface  macrogeometry  coefficients 

53.  The  surface  macrogeometry  coefficients  S for  wedges^  sire 

plotted  in  plate  43  as  20  log  S versus  wedge  interior  angle  0 ; 20 

log  S for  steps  are  plotted  in  plate  44  versus  2nfH/V  , where  f 

s 

* frequency,  Hz;  H = step  height,  m;  and  V = shear  wave  velocity,  m/sec. 

s 

A wedge  is  defined  by  two  plane  surfaces  that  Join  to  form  a change  in 
slope.  A step  is  an  abrupt  change  in  elevation  and  is  really  formed  by 
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two  wedges  at  90  deg,  separated  in  distance  by  the  height  of  the  step. 
Normally,  surface  macrogeometry  coefficients  can  be  assumed  to  be  unity 
(20  log  S = 0)  unless  the  terrain  feature  is  large  in  comparison  with  a 
wavelength  of  the  propagating  wave,  i.e.  the  feature  is  larger  than, 
for  example,  a half  wavelength,  X (in  meters)  where  X = V^/f  . 

54.  The  wavelengths  of  propagating  waves  are  a function  of  fre- 
quency and,  for  wide-band  sources,  would  have  wide  ranges  of  values. 

Thus,  it  might  be  advisable  to  make  a prediction  for  a site  as  though  it 
were  flat  terrain  to  determine  the  site  predominant  frequency.  If  the 
wavelength  of  the  predominant  frequency  were  large  compared  with  the 
macrogeometry  feature , the  effect  would  be  relatively  small;  however, 
even  small  terrain  features  will  have  an  effect  on  the  high-frequency 
components.  Thus,  if  accurate  predictions  are  desired,  complex  terrain 
features  must  be  described  by  using  a combination  of  wedges  and  steps 
and  summing  the  effects  of  each.  For  steps,  the  value  of  20  log  S can 
be  determined  (see  plate  44)  as  a function  of  step  height;  for  wedges, 
the  distance  over  which  the  wedge  is  formed  must  be  large  with  respect 
to  a wavelength  before  the  value  of  20  log  S need  be  applied,  and  the 
appropriate  value  can  be  selected  from  plate  43. 

55-  The  macrogeometry  coefficients  for  a wedge  are  functions  of 
the  wedge  angle  6 (measured  as  an  angular  deviation  from  the  approach- 
ing plane  surface)  and  Poisson's  ratio  a . Poisson's  ratio  can  be 
defined  in  terms  of  compression  and  shear  wave  velocities,  and  Vg  , 

as 


1 - 2<  VV2 

2 - 2 (v./vj2 

S p 


(10) 


Waves  propagated  in  either  direction  on  a wedge  have  the  same  macro- 
geometry coefficients. 

56.  In  contrast  to  the  wedge  coefficients,  the  step  coefficients 
are  functions  of  frequency;  as  frequency  increases,  the  amplitude  of 
the  wave  traversing  a step  is  reduced  by  the  coefficient.  Also,  two 
curves  are  shown  for  the  step  coefficients;  one  for  the  wave  incident 
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to  the  elevated  side,  and  one  for  the  wave  incident  on  the  depressed 
side. 

57*  The  macrogeometry  coefficients  are  given  for  homogeneous 
half-space  materials  only.  However,  as  a first  approximation,  the  homo- 
geneous half-space  coefficients  can  he  used  for  a layered  site,  using 
the  material  properties  of  the  layer  that  carries  most  of  the  seismic 
energy. 


Graphic  Modeling  Procedure 

58.  Three  examples  are  given  to  illustrate  the  graphic  modeling 
procedure.  This  procedure  requires  the  user  to  select  one  curve  each 
for  the  source  stress  signal,  the  source  coupling  coefficients,  the 
transmission  coefficients,  and  the  macrogeometry  coefficients.  Values 
from  each  of  the  curves  are  selected  at  various  frequencies  and  tabu- 
lated. The  microseismic  signal  is  then  constructed  by  adding  the  values 
at  each  frequency  and  plotting  the  results  versus  frequency. 

Site  1 results 

59*  Particle  velocity  signals  were  predicted  for  a man  walking 
at  a range  of  10  m on  a 0.75-m-thick  layer  of  material  B over  an 
infinite  half-space  of  material  F.  A source  stress  signal  representing 
a footstep  on  a hard  surface  under  a low  compaction  condition  (plate  2) 
was  selected  to  demonstrate  signal  propagation  over  a flat  surface 
with  approximately  the  same  properties  as  the  subsurface  material. 

Table  2 summarizes  the  values  read  from  the  figures  and  plates  identi- 
fied at  the  bottom  of  the  table  at  the  frequencies  indicated  in  the 
left-hand  column.  The  graphic  solution  is  shown  in  fig.  6 together 
with  a digital  computer  solution  of  the  same  problem.  Agreement  be- 
tween the  two  solutions  is  excellent  at  most  frequencies,  except  near 
the  lower  bound  where  the  digital  computer  solution  gives  somewhat 
larger  values  than  the  graphic  solution.  The  digital  computer  time- 
domain  prediction  for  particle  velocity  shown  in  fig.  7 is  the  inverse 
Fourier  transform  of  the  digital  computer  frequency-domain  signal  in 
fig.  6.  This  time-domain  signal  is  shown  here  as  a matter  of  interest 
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lple  site  1:  frequency-domain  signals  over  flat  surface  of  a 0.75-m-thick  layer  of 

. B on  an  infinite  half-space  of  material  F;  source  stress  signal  from  plate  2 


inrple  site  1:  time-domain  signal  over  flat  surface  of  a 0.75-^-thick  layer  of  material 
over  an  infinite  half- space  of  material  F;  source  stress  signal  from  plate  2 


(as  are  the  rest  of  the  time-domain  signals  shown  in  this  section), 
since  the  graphic  pr  .jess  results  in  a frequency-domain  signal  only. 

Site  2 results 

60.  Particle  velocity  signals  were  predicted  for  a 3-m  layer  of 
material  B over  an  infinite  half-space  F in  the  same  manner  as  shown  for 
site  1 in  table  2.  A source  stress  signal  representing  a footstep  on 

a flat  hard  surface  under  a low  compaction  condition  (plate  2)  was  se- 
lected as  before.  Note  that  except  for  the  layer  thickness  (i.e.  3m 
instead  of  0.75  m) , the  conditions  are  the  same  as  for  the  conditions 
at  site  1.  The  graphic  solution  for  ranges  of  5,  10,  20,  30,  and  40  m 
is  shown  in  fig.  8.  The  digital  computer  time-domain  predictions  for 
particle  velocity  corresponding  to  the  graphic  predictions  in  the  fre- 
quency domain  are  shown  in  fig.  9 • When  the  10-m-runge  graphic  solu- 
tions of  figs.  6 and  8 are  compared,  a broad  peak  with  a predominant 
frequency  of  approximately  50  Hz  is  shown  for  the  3-m-thi ck-layer  site; 
whereas  the  0.75-m-thick-layer  site  has  a slightly  narrower  peak  with  a 
predominant  frequency  of  80  Hz.  In  the  time  domain  the  3-m-thick-layer 
site  (fig.  9)  has  a high-amplitude  short-duration  signal,  and  the  0.75-*n- 
thick-layer  site  (fig.  7)  has  a lower-amplitude  long-duration  signal. 
Also,  as  seen  in  fig.  8,  as  the  range  increases  (i.e.  from  5 to  40  m) , 
the  predominant  frequency  of  the  site  shifts  from  60  to  25  Hz,  and  the 
amplitude  at  the  predominant  frequency  decreases  with  range. 

Site  3 results 

61.  Site  3 is  identical  with  site  2 in  substrate  conditions,  i.e. 
they  are  both  composed  of  a 3.0-m  layer  of  material  B over  infinite  half- 
space of  material  F.  However,  the  surface  material  at  site  3 is  not  as 
firm  as  that  at  site  2.  The  source  stress  signal  for  this  "medium" 
surface  under  plastic  conditions  was  selected  from  plate  3 (case  3). 

Also,  site  3 exhibits  an  abrupt  change  in  slope;  whereas  site  2 is  flat. 
This  change  in  slope  is  represented  by  a 50-deg  wedge,  and  its  effect  on 
a propagating  signal  can  be  obtained  from  plate  U3.  In  this  example,  a 
Poisson's  ratio  of  0.4  was  used;  therefore,  the  value  20  log  S is  6. 

62.  Fig.  10  shows  the  computer-generated  relations  of  particle 
velocity  versus  time  for  site  3 at  a range  of  10  m.  These  relations  do 
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so 


Example  site  2:  frequency-domain  signals  over  flat  surface  at  various  ranges,  3.0-m- 

• of  material  B on  infinite  half-space  of  material  F;  source  stress  signal  from  plate 


material  B on  infinite  half-space  of  material  F;  source  stress  signal  from  plate 


not  account  for  the  effects  of  the  surface  macrogeometry  feature,  and, 
therefore,  can  be  compared  directly  with  those  in  plot  b,  fig.  9,  to 
estimate  the  effect  of  changing  the  input  stress  signal  (i.e.  change 
from  case  2 to  case  3,  see  plates  2 and  3)  on  the  particle  velocity-time 
relations.  The  softer  surface  material  resulted  in  a reduction  in  peak 
amplitude  from  5 to  2.5  * 10“^  cm/sec  and  a signal  of  noticeably  longer 
duration.  Also,  the  predominant  frequency  for  site  3 was  less  than  that 
predicted  for  site  2 (i.e.  about  20  compared  with  50,  see  figs.  8 and 
10).  The  effect  of  the  signal  propagating  over  the  50-deg  wedge  can 
also  be  seen  in  fig.  10.  The  result  was  obtained  by  subtracting  6 db 
from  the  prediction  of  propagation  over  the  flat  surface;  6 db  is  equi- 
valent to  a 50  percent  reduction  in  the  amplitude  of  the  particle 
velocity-time  relation  shown  in  fig.  11. 

Summary  and  Recommendations 

63.  The  study  has  primarily  dealt  with  the  development  and  pre- 
sentation of  simple  graphic  procedures  that  can  be  used  to  predict 
seismic  signal  amplitudes  and  frequency  characteristics.  During  this 
phase  of  the  study,  exhaustive  analysis  of  the  relations  derived  to 
develop  the  graphic  procedures  was  not  accomplished;  nevertheless,  the 
following  generalizations  can  be  advanced. 

Summary 

6U . A simple  but  useful  graphic  solution  to  the  wave  generation 
and  propagation  phenomena  for  man-walking  targets  was  developed.  The 
solution  is  such  that  the  terrain  effects  on  (a)  target-ground  interac- 
tion, (b)  energy  source  coupling,  (c)  transmission,  and  (d)  surface 
macrogeometry  can  be  isolated  and  studied,  thereby  permitting  the  user 
to  determine  quantitatively  how  the  separate  contributors  (a-d  above) 
interact  synergistically.  Generalizations  based  on  the  curves  generated 
in  this  report  must  be  made  within  the  parameter  ranges  used  to  generate 
the  curves.  It  is  emphasized  that  whenever  the  generalizations  are  ap- 
plied in  situations  outside  the  parameter  ranges  from  which  extractions 
were  made,  errors  may  arise.  However,  interpolations  within  the 
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.1.  Example  site  3:  time-domain  signals  over  flat  and  sloped  surfaces,  3-m- thick  layer  of 

material  B over  infinite  half-space  of  material  F;  source  stress  signal  from  plate  3 


parameter  ranges  studied  can  be  made  with  confidence,  and,  therefore, 
generalizations  concerning  cases  not  actually  studied  in  this  report  are 
possible  (paragraphs  9 and  10). 

65.  To  mathematically  simulate  the  wave  generation  and  propagation 
phenomena,  a soil  model  consisting  of  at  least  the  following  factors 

is  required  (paragraphs  5-8): 

a.  Force-deflection  (soil  surface)  relations. 

b.  Seismic  properties  of  each  propagating  medium  (Vs  , Vp  , 

Y , and  layer  thickness)  to  a depth  approximately  equal 
to  a half  wavelength  of  the  lowest  frequency  of 
importance. 

£.  The  range  from  signal  source  to  sensor. 

d_.  Generalized  surface  macrogeometry  features,  i.e.  fea- 
tures described  as  combinations  of  steps  and  wedges. 

_e.  Soil  damping  factor  (in  this  study  an  average  damping 
factor  of  0.03  was  used). 

66.  The  results  of  this  study  show  (paragraphs  40-43*  48-52,  and 
58-62)  that  site  properties  influence,  to  a large  degree,  the  character 
of  the  signal  from  a footstep.  The  signal  arriving  at  the  sensor  is 
affected  by  (a)  target-ground  interaction,  (b)  energy  coupling  to  the 
substrate,  (c)  propagation  from  the  source  to  the  sensor,  and  (d)  trans- 
mission of  the  wave  over  surface  geometry  features.  The  graphic  solu- 
tion developed  as  depicted  in  fig.  4 requires  the  summation  of  the  ef- 
fects enumerated  above.  Because  of  this  summation  process  over  the 
frequency  range  of  interest  it  has  been  established  that,  if  the  source 
signal  has  certain  frequency  components  that  have  extremely  low  ampli- 
tudes, these  components  will  have  large  negative  logarithmic  values  and 
will  not  contribute  substantially  to  the  total  seismic  signal.  The  only 
way  these  frequency  components  can  form  a part  of  the  signal  is  by  way 
of  the  uncontrollable  background  noise.  In  the  same  manner,  if  the 
source  coupling  or  the  site  transfer  function  has  low  amplitude  over 
parts  of  the  frequency  spectrum,  these  parts  will  be  highly  attenuated 
in  the  seismic  signal  that  is  generated. 

Recommendations 

67.  Inasmuch  as  the  analytical  tools  for  studying  the  effects  of 
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terrain  on  Rayleigh  wave  generation  and  propagation  are  Just  emerging 
from  research,  experience  gained  from  use  of  the  models  is  still  quite 
limited.  For  this  reason,  generalizations  or  "theorems"  useful  to  the 
design  and  deployment  of  UGS  must  be  derived  through  deliberate  applica- 
tion of  the  models.  Therefore,  it  is  recommended  that  results  of  the 
study  reported  herein  be  considered  interim,  and  that  the  study  of  the 
terrain  effects  on  generation  and  propagation  of  seismic  signals  from 
man-walking  targets  be  continued.  It  is  especially  desirable  that  a 
more  comprehensive  set  of  signal-terrain  relations  be  generated  and 
studied  with  the  objective  of  collapsing  the  myriad  of  curves  to  a 
simpler  and  more  widely  applicable  graphic  solution. 

68.  It  is  further  recommended  that  the  study  be  expanded  to  in- 
clude multiple  personnel  and  vehicle  targets.  It  is  particularly  de- 
sirable that  forcing  functions  (stress-time  signals)  be  derived  for  a 
family  of  Warsaw  Pact  vehicles,  and' that  these  forcing  functions  be  used 
in  the  sensitivity  studies.  It  is  believed  that  the  work  recommended 
above,  in  conjunction  with  data  presented  in  this  report,  will  be  ex- 
ceptionally useful  to: 

a.  Fill  in  gaps  in  sensor  design  data  base. 

b.  Identify  target  feature  consistency  (i.e.  features 
insensitive  to  terrain  variations ) . 

c_.  Design  sensor  fields  (i.e.  calculate  effective  range 
of  devices  to  determine  layout  and  spacing  for  optimum 
protection ) . 

d.  Simplify  site  requirements  for  mapping  criteria  (i.e. 
the  model  can  be  used  to  identify  site  parameters  that 
' can  be  combined  or  inferred  from  other  more  easily 
identifiable  site  characteristics). 
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Table  1 

Seismic  Properties  of  Media 


Phase  Velocities,  m/sec 


Material 

Shear 

Wave 

V 

s 

Compression 

Wave 

V 

P 

Rayleigh 

Wave* 

VR 

Density 

0 

g/cm 

A 

40 

ICO 

38 

1.3 

B 

120 

300 

113 

1.6 

C 

200 

500 

189 

1.9 

D 

200 

500 

189 

2.0 

E 

500 

1000 

467 

2.3 

F 

1000 

2000 

937 

2.6 

G 

4oo 

1000 

377 

1.55 

* Rayleigh  wave  phase  velocity  refers  to  that  for  homogeneous  half- 
space only;  for  layered  half -space,  refer  to  curves  of  Vr  versus 
frequency. 
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In  accordance  with  ER  70-2-3,  paragraph  6c(l)(b), 
dated  15  February  1973*  a facsimile  catalog  card 
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